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Abstract. The explanation of the suppression of high-pT hadron yields at RHIC in terms of jet-quenching 
implies that the multiplicity distributions of particles inside a jet and jet-like particle correlations differ 
strongly in nucleus-nucleus collisions at RHIC or at the LHC from those observed at e + e~ or hadron 
colliders. We present a framework for describing the medium-induced modification, which has a direct 
interpretation in terms of a probabilistic medium-modified parton cascade, and which treats leading and 
subleading partons on an equal footing. We show that our approach can account for the strong suppres- 
sion of single inclusive hadron spectra measured in Au-Au collisions at RHIC, and that this implies a 
characteristic distortion of the single inclusive distribution of soft partons inside the jet. We determine, as 
a function of the jet energy, to what extent the soft fragments within a jet can be measured above some 
momentum cut. 

PACS. 12.38.Mh Quark-gluon plasma - 25.75.-q Relativistic heavy-ion collisions 
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1 Introduction 

Among the most notable results from the first years of 
nucleus-nucleus collisions at RHIC stand the deficit in 
high transverse-momentum hadrons and the suppression 
of leading back-to-back hadron correlations observed in 
central Au-Au collisions with respect to expectations from 
scaling the yields measured in pp or d-Au collisions pQ. 
These observations are consistent with the "jet-quenching" 
scenario: before they hadronize in the vacuum, partons 
produced in the dense matter created in head-on Au-Au 
collisions dissipate a significant fraction of their energy 
through an enhanced radiation of soft gluons ■ 

Irrespective of the details of the implementation of the 
medium-enhanced radiation of gluons — either through 
coherent multiple soft-momentum transfers a la Landau- 
Pomcranchuk-Migdal 5,6 , or through single hard scat- 
tering |3] — jet-quenching models of inelastic (radiative) 
energy loss are quite successful in explaining present light- 
hadron data from RHIC [312111] ■ Nevertheless, there 1 re- 
mains much room for technical improvement over the ex- 
isting formulations of inelastic energy loss. Thus, a generic 
feature of these approaches is that they only consider 
the medium-induced enhancement in gluon radiation for 
the leading parton, discarding the medium influence on 
subleading partons. Such an approximation may remain 
under control when dealing with leading-hadron produc- 
tion; however, predictions involving subleading particles 
become questionable, be it for jet shapes, which may be- 
come experimentally accessible at the LHC ^U] , or simply 
for intrajet two-particle correlations. Similarly, in exist- 



ing models energy-momentum conservation is not explic- 
itly conserved at each parton splitting, but only globally, 
through various ad hoc corrections. 

A novel formulation of medium-induced parton en- 
ergy loss was recently introduced in Ref. which aims 
at correcting some of the shortcomings of standard ap- 
proaches mentioned above. Thus, it is the first one that 
deals equally with the various splittings of both leading 
and subleading partons inside a shower. Furthermore, it 
automatically conserves energy-momentum at each parton 
splitting. These improvements are admittedly obtained, 
at the moment, at the cost of a less accurate treatment 
of other aspects of radiative energy loss, which should be 
restored in a future Monte-Carlo implementation. Never- 
theless, the physical constraints we choose to emphasize 
in our approach may prove to be of greater importance 
than those retained in other existing models, in particu- 
lar for the discussion of high-px particle correlations (see, 
however, Refs. 1 1 211 1 3| for an alternative extension of par- 
ton energy loss which allows one to tackle that issue) or 
multiplicity distributions inside jets. 



2 Formalism 

One of the most testing ground of the color structure of 
quantum chromodynamics (QCD) is provided by the jets 
that are created in elementary-particle collisions, be it in 
e + e~ or in pp/pp collisions. Thus, the asymptotic shape 
of the distribution of hadron momenta inside a jet can 
be computed exactly, especially at small momentum frac- 
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tions x = p/Ejct, by a resummation of infrared-singular 
terms to all orders, within the so-called Modified Lead- 
ing Logarithmic Approximation (MLLA) of QCD [1411151 
ITT)] . Color coherence thus results in destructive interfer- 
ence between partons, which leads to a suppression of 
hadrons with small x. The interference is actually equiva- 
lent, to double and single logarithmic accuracy in ln(l/cc) 
and hx{Q/A e s) — where Q ~ Ej Ct is the jet virtuality 
and A c ff an infrared cutoff which is eventually fitted to 
experimental data — to an angular ordering of the se- 
quential parton decays within the shower, with leading- 
order splitting functions. An important prediction of this 
angular-ordered probabilistic parton cascade is, to next- 
to-leading order <Jas, the characteristic "hump-backed 
plateau" shape of the distribution of parton momenta in- 
side a jet, represented as a function of ln(l/x). The parton 
shower, evolved down to an infrared cutoff Qo ~ A e g, is 
eventually identified to a hadron jet, by mapping locally 
each parton onto a hadron ( "Local Parton-Hadron Dual- 
ity", LPHD): for each hadron type, the hadron distribu- 
tion equals K h times the parton distribution, where K h is 
a proportionality factor of order unity. 1 This resummation 
and the LPHD prescription give a good description of the 
measured longitudinal distributions of hadrons D h (x, Q 2 ) 
over a wide energy range, both in e + e~ |17I18| and in 
PP PHI collisions. For instance, Fig. ^ shows the distribu- 
tion of inclusive hadrons D h (x,Q 2 ) inside 17.5 GeV jets, 
a typical jet energy at RHIC, as measured in e + e~ anni- 
hilation by the TASSO Collaboration [Uj, together with 
the MLLA prediction with K h = 1.35. ' 

The formalism developed in Ref. ^JJ to describe the 
medium-induced distortion of jets reduces to the MLLA 
baseline in the absence of a medium. This new approach 
involves different approximations, emphasizing other phys- 
ical aspects, from the standard models of parton energy 
loss that are currently used in the phenomenology of RHIC 
data. Thus, present model comparisons to RHIC data 
start with a medium-modified energy spectrum of radi- 
ated gluons, d/ tot = d/ vac + d/ med pSEj . The part corre- 
sponding to the "normal" vacuum radiation shows a dou- 
ble logarithmic dependence d/ vac = ff^r^r; its integral 
over k gives rise to the leading InQ 2 term in the DGLAP 
evolution equation. This contrasts to the k-integration 
of d/ med , which is infrared- and ultraviolet-safe |H| and 
leads to a nuclear-enhanced "higher-twist" contribution, 
oc qL/Q 2 , where q is the transport coefficient that charac- 
terizes the medium, subleading in an expansion in l/Q 2 , 
but enhanced with respect to other such terms by a fac- 
tor proportional to the geometrical extension ~ L of the 



1 In addition to using different phenomenological proportion- 
ality constants K h , the hadron- type dependence in MLLA is 
also implemented by cutting the development of the parton 
cascade at different scales Qo, with Qo chosen equal to the 
hadron mass. This requires using a slightly more elaborate form 
of the spectrum than Eq. Q, see Ref. |16| . and leads to an 
energy-independent stiffening of spectra for heavier hadrons. 
The comparison to inclusive yields in elementary-particle col- 
lisions were done with a unique cutoff parameter for all particle 
types, Qo = A cS = 253 MeV. 




£ = ln(l/jc) 



Fig. 1. Longitudinal distribution diV/dln(l/a;) of inclusive 
hadrons inside a jet of energy Ej ot = 17.5 GeV, as a function of 
ln(l/x) = ]n(Ej e t/p), as measured by TASSO ^7] and within 
MLLA (solid curve: / me d = 0; dashed curve: / mc d = 0.8). 

target. In practice, however, the parton virtuality does 
not enter the existing comparisons to experimental data, 
where one rather considers the k-integrated gluon distri- 
bution u) d ^2° , neglecting the Q 2 -dependence. In addi- 
tion, existing approximations only include the extra source 
of gluon radiation d/ med for the leading parton, dropping 
it for the further medium-induced splittings of subleading 
partons in the shower. 

The obvious way to improve over this state of the art 
is to replace the double differential gluon spectrum d/ vac 
by d/ tot in all leading and subleading splitting processes 
of a medium-modified parton cascade. This can only be 
done within a Monte-Carlo approach, which we intend 
to develop in future studies. The first step in that direc- 
tion, which is still fully analytical, consists in using an 
additional approximation: instead of using the computed 
k-intcgrated medium-induced distribution ui ^ — , we re- 
placed it by a constant / m ed- One can check that this 
is not too crude an assumption in the kinematic regime 
tested at RHIC, where it amounts to a similar uncertainty 
as that arising from whether one should use the multiple- 
soft scattering approach or the single hard scattering pic- 
ture. Under this assumption, we have then used the extra 
medium-induced spectrum ui d ^ on the same level as 
to d 4— , i.e., as a leading logarithmic correction JJJ. With 
this ansatz, our formalism ensures energy- momentum con- 
servation at each parton splitting, and treats all leading 
and subleading parton splittings on the same footing. 



3 Phenomenological predictions 

As stated above, in our analytical approach we approxi- 
mate the medium-induced contribution to the gluon spec- 
trum u> dI di ^ by a constant / mo d in the kinematically rel- 
evant range of to. This actually amounts to considering 
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"medium-modified parton splitting functions" , which dif- 
fer from the standard ones by enhancing their singular 
parts by a factor (1 + / me d)- 2 This formulation allows us 
to follow the same line of technical arguments as that used 
for the calculation of jet multiplicity distributions in the 
absence of a medium [TB] , 



3.1 Distortion of the hump-backed plateau 

We can thus compute the momentum distribution of par- 
tons within a (gluon-initiated) parton cascade: 



xD g (x,r)- 



4JV C (1 + / med )r 



bB(B + l) 
r+io ° dv 

I — : X~ v ${-A + B + 1,5 + 2;-vt), 



2vri 



(1) 



where r = 
function, 
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ln(Q/A c ff), <P is the confluent hypergeometric 



A 

1^ /med 



4iV c (l + / med ) 



bv 

Nc+ 3W 



R - a 

h = l ±N c - 2 -N 
3 3 



f> 



and Nf is the number of active flavors. Quite naturally, 
by setting / mo d = in Eq. QJ one recovers the usual 
"vacuum" expression of x D g {x, r). 

To exemplify the effect of the medium-enhanced gluon 
radiation on the hump-backed plateau of particle produc- 
tion, we compare in Fig.^the spectra Eq. JJJ for a jet with 
energy E ict — 17.5 GeV, shown as a function of \n{l/x), 
in the cases / mo d = (no medium) and / me d = 0.8, where 
the choice of the latter value will be justified in Sec. 13. 31 3 
One clearly sees that the effect of the medium is a strong 
distortion of the distribution, with a depletion of the num- 
ber of particles at large x, and correspondingly a largely 
enhanced emission of particles at small x: due to energy- 
momentum conservation in the parton cascade, the energy 
which in the vacuum is taken by a single large- x parton is 
redistributed over many small- x partons in the presence 
of a medium. 

An alternative way to illustrate the influence of the 
medium on the spectra is to display the distributions as a 
function of the momentum fraction x, see Fig.El The plot 
emphasizes the suppression in high- a; hadron production, 
which is the aspect tested by experimental studies of the 
high-px spectra of hadrons, rather than the enhancement 
in small- a; particles. 

2 Such a modification of parton splitting functions was dis- 
cussed in Ref. 1201 . where it results from considering nuclear- 
enhanced twist-four parton matrix elements in studies of 
deeply inelastic eA scattering. 

3 Both spectra were actually scaled by a factor K h — 1.35, so 
that the vacuum one match the corresponding TASSO data. 
This goes along the idea that the medium affects the devel- 
opment of the parton shower, while hadronization takes place 
"in the vacuum" and thus the LPHD coefficient should be the 
same in both cases. 
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Fig. 2. Single inclusive distribution of hadrons, dN/dx, inside 
a jet of energy E- 3e % = 17.5 GeV, as a function of x = -Ejet, 
as measured by TASSO |17) and within MLLA (solid curve: 
/med = 0; dashed curve: / mcd = 0.8). 



3.2 Associated multiplicity 

Once the longitudinal multiplicity distribution inside a jet 
is known, a straightforward integration yields the number 
of hadrons inside the jet with transverse momenta larger 
than a given cut, N(pr > P™')- We- can then calculate 
this multiplicity for jets with the same energy both in 
the presence of medium effects (in which case, the lower 
cut gives some control on the high-multiplicity soft back- 
ground of nucleus-nucleus collisions at RHIC or LHC, 
over which the jet develops) and in vacuum, and compute 
their ratio. We display the latter in Fig. [3] for jets with 
-Ejet = 17.5 GeV and medium effects modeled by a con- 
stant coefficient / me d = 0.8. The ratio is seen to be smaller 
than 1 for p™* > 1.5 GeV/c, while the medium-induced 
enhancement in soft-particle production becomes domi- 
nant for smaller values of the transverse- momentum cut. 
The crossover value is close to that reported by the STAR 
Collaboration in attempts at measuring the excess of par- 
ticles inside the back jet over the soft background \21\. 

A more thorough comparison between our calculation 
and the STAR data would have to take into account sev- 
eral effects which we did not consider in our simple ap- 
proach. First, in our calculation the energy of the jet 
is known and fixed, while experimentally it spans some 
range in energy. In the same line, we did not try to re- 
produce the experimental fact that the high-p^ particles 
measured at RHIC come from jets that have traversed 
different in-medium path lengths, and thus radiated dif- 
ferent energy fractions. Eventually, we did not model the 
soft background of Au-Au collisions at RHIC, which ex- 
tends well over the region where the medium-enhanced 
part of the multiplicity lies. Nevertheless, the reasonable 
agreement between our calculation and the experimental 
measurement is a further hint that energy is indeed redis- 
tributed from high- to low- a; partons through the influence 
of the medium, although nothing definitive on the actual 
shape of the distorted spectrum of radiated partons can 



4 



Nicolas Borghini: Multiplicity distributions inside parton cascades developing in a medium 




Fig. 3. Ratio of the multiplicities of particles above a lower 
momentum cut inside a jets of energy s£j e t = 17.5 GeV 
developing in a medium (/ me d = 0.8) and in the vacuum. 



be claimed at the moment. Note, however, that the cal- 
culations reported in Ref. indicate that the crossover 
between enhancement and depletion should take place at 
transverse momenta p™' ~ 4 — 7 GeV/c for jets of energy 
-Ejet = 100 — 200 GeV, which should be accessible at the 
LHC. This should leave a window above the upper kine- 
matic boundary of the soft background, in which there is 
an enhancement of the jet multiplicity, thereby allowing 
a more detailed characterization of the medium-enhanced 
radiation. 



3.3 Nuclear modification factor 

As a last phenomenological prediction of our model, let 
us compute the medium- modification factor Raa, defined 
as the ratio of the hadron yield in nucleus-nucleus over 
that in pp collisions, scaled by the number of binary colli- 
sions. For that purpose, we need to supplement a further 
ingredient to our approach, namely the underlying par- 
tonic momentum distribution at RHIC energies. Since we 
only want to show that our implementation of radiative 
parton energy loss is able to result in a factor of 5 suppres- 
sion of the hadron yields, as measured experimentally, we 
shall restrict ourselves to using a power-law parton spec- 
trum We do not take into account the kinematic- 
boundary influence on the shape of the pT-spectrum, nor 
"initial-state" effects like shadowing or Cronin effect. Even- 
tually, as in Sec. 13.21 we do not attempt to model the ge- 
ometry of the medium, nor the fluctuations in the energy 
radiated by partons with a given energy. 

The result of our calculation is shown in Fig. 0] to- 
gether with the PHENIX ir° data at = 200 GeV 
We see that our choice of / me d = 0.8 for the coefficient that 
models medium effects allows us to reproduce a modifica- 
tion factor of ~ 0.2, which justifies our using it in all 
applications reported here. This suppression by a factor 5 
actually reflects the observation one can make on Fig. [21 
that the spectrum of radiated partons in a medium is a 
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Fig. 4. Transverse-momentum dependence of the nuclear mod- 
ification factor Raa, computed for a medium-enhanced radia- 
tion of gluons modeled by / mc d = 0.8, compared to PHENIX 
data on tt° production [22| 

factor ~ 3 — 6 times smaller than that of partons within a 
jet in vacuum for momentum fractions x ~ 0.4—0.6, where 
Raa mostly tests the spectra. Such values of x are often 
considered to lie outside the small- a; region where MLLA 
is a proper approach to hadron production in elementary- 
particle collisions. Nevertheless, the possible uncertainty 
of using such a resummation of large ln(l/x) terms (in- 
stead of, for instance, fitted fragmentation functions) in 
computing the nuclear modification factor is relatively 
small, inasmuch as Raa is a ratio, which does not depend 
on the overall normalization of the vacuum spectrum, and 
that we are discussing a suppression by a factor of 5, not 
an effect of a few percent. 

4 Conclusion 

We have reported a first step towards a description of 
parton cascades developing in a medium, which conserves 
energy-momentum at each successive parton splitting, and 
treats all partons in the shower on the same footing 
The simplified analytical formalism we have presented is 
able to reproduce at least semi-quantitatively several char- 
acteristic features of RHIC data, such as the suppression 
of high-momentum particle yields and the enhanced soft- 
particle distribution associated to high-p^ trigger parti- 
cles. The price we pay for an analytical treatment is the 
inaccurate handling of the Q 2 ~dependence of parton ra- 
diation, which is not dealt with more properly in other 
existing formalisms of radiative parton energy loss, and 
the approximate treatment of the w-dependence of the 
spectrum, which might not be of prime importance in the 
RHIC regime. 

In the future, we expect to remove the latter deficien- 
cies by implementing our approach within a Monte-Carlo 
formulation. In that respect, the analytical results pre- 
sented here will serve as a powerful consistency check for 
the correct implementation of splitting processes in the 
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Monte-Carlo simulation. Such simulations will allow us to 
address the issue whether and how medium-induced par- 
ton energy loss depends on virtuality, as heuristic argu- 
ments suggest it should: since a hard parton of virtuality 
Q and energy E has a lifetime ~ E/Q 2 in the rest frame 
of the dense matter through which it propagates, before 
it branches, high- virtuality partons will branch so quickly 
that their splitting is not influenced by the presence of 
the medium (so that / mc d should decrease with increasing 
Q 2 ). This question will be particularly important for ad- 
dressing the logarithmically wide kinematic range acces- 
sible in Pb-Pb collisions at the LHC. An angular-ordered 
Monte-Carlo parton cascade will also shed light on other 
issues that have been studied in the absence of medium 
effects, as two-particle correlations within jets 115023124) . 
transverse jet broadening, or the differences between light- 
and heavy-quark jets [251 - Finally, it will also allow one to 
implement realistic features of the collision, as e.g. fluc- 
tuations or the geometry of the dense created medium 
traversed by the fast parton. 
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